Effects of spin-orbit interaction on superconductor-ferromagnet heterostructures: 
spontaneous electric and spin surface currents 
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We find proximity-induced spontaneous spin and electric surface currents, at all temperatures 
below the superconducting T c , in an isotropic s-wave superconductor deposited with a thin ferro- 
magnetic metal layer with spin-orbit interaction. The currents are carried by Andreev surface states 
and generated as a joint effect of the spin-orbit interaction and the exchange field. The background 
spin current arises in the thin layer due to different local spin polarizations of electrons and holes, 
which have almost opposite velocities in each of the surface states. The spontaneous surface electric 
current in the superconductor originates in asymmetry of Andreev states with respect to sign rever- 
sal of the momentum component parallel to the surface. Conditions for electric and spin currents 
to show up in the system, significantly differ from each other. 
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PACS numbers: 74.45.+C, 74.81.-g 

Proximity effects in superconductor-ferromagnet het- 
erostructures have attracted much attention for recent 
years. In contrast with the nonmagnetic case, magnetic 
surfaces and interfaces make spin-flip processes possible 
and suppress an s-wave superconducting order param- 
eter, generating Andreev bound states in adjacent su- 
perconducting region o 1 ! 2 ! 3 ! 4 ! 5 . Spin structure of Andreev 
bound states near complex magnetic interfaces can be 
rather involved 5 . Triplet components of the order param- 
eter in a singlet superconductor can be induced by fer- 
romagnets under certain conditions 6 :. Cooper pair wave 
functions exponentially decay into the bulk of ferromag- 
nets, oscillating at the same time 7 -, and acquire a triplet 
component in the ferromagnetic region 8 . These prox- 
imity effects can lead, in particular, to specific proper- 
ties of the Josephson current through magnetic inter- 
faces, which have been intensively studied both theoret- 
ically and experimentally 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 7 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 . Also, 
proximity-induced nonmonotonic dependence of the su- 
perconducting critical temperature on the thickness of 
the ferromagnetic layer has been thoroughly studied 
for superconductor-ferromagnetic metal bilayers or het- 
erostructures (see, for example ^ILi^i^iLl^ and refer- 
ences therein). 

Spontaneous surface currents represent other impor- 
tant example of possible proximity-induced effects. Spon- 
taneous electric currents, taking place near surfaces or 
interfaces on the scale of the superconducting coherence 
length, produce a magnetic field and, hence, a counter- 
flow of screening supercurrents on the scale of the pen- 
etration depth. The electric current can arise, for ex- 
ample, near nonmagnetic surfaces and interfaces of un- 
conventional superconductors, whose states break time- 
reversal symmetryiSiSS. In particular, the electric cur- 
rent, carried by Andreev states, appears near nonmag- 
netic surfaces and interfaces of chiral superconductors^. 
The current also occurs, if a surface-induced subdomi- 
nating pairing shows up near surfaces of d-wave super- 
conductors, breaking time-reversal symmetry of the su- 
perconducting state 22 . Other possible mechanism gener- 



ating electric surface currents, is specifically based on a 
paramagnetic response of the zero-energy Andreev sur- 
face states to an applied magnetic field. This can take 
place at low temperatures at smooth (110) surfaces of 
d-w&ve superconductors 2 ^ 2 ^, as well as in a system with 
a thin ferromagnetic metallic layer deposited on a semi- 
infinite bulk isotropic s-wave superconductor—. In the 
latter case the energy of the surface states becomes zero 
only for several values of the layer thickness and, in 
the presence of particle-hole asymmetry, the spontaneous 
electric current is accompanied by a spontaneous surface 
spin current. Dissipationless background spin currents, 
which take place in various systems in the equilibrium 
and do not lead to any spin accumulation, have been 
a subject of recent discussions and studies2SiSL22. The 
spin currents can be generated, for example, by the spin- 
orbit interaction (in particular, via Rashba term) in two- 
dimensional metals. Measurements of background persis- 
tent spin currents are not carried out for now, although 
some suggestions for a direct detection of these currents 
have been proposed in the literatur e 26 ' 29 . 

In the present paper we study spontaneous currents 
under conditions, when spin-orbit interaction takes place 
in a thin ferromagnetic metal layer in proximity to an 
isotropic s-wave superconductor. Joint effect of the spin- 
orbit interaction, described by the Rashba term, and the 
exchange field is shown to play an important role in gen- 
erating spontaneous currents. We find that superconduc- 
tor induces background spin currents in the ferromag- 
netic layer with the spin-orbit interaction (FSOL) for all 
temperatures below the superconducting T c . This spin 
current is carried by Andreev surface states and takes fi- 
nite values due to different local spin polarizations of elec- 
trons and holes, which have almost opposite velocities in 
each of these states. Maximal possible values of the back- 
ground spin current density are the order of the Landau 
depairing current density. We find proximity-induced 
finite spin currents within the quasiclassical approach, 
when only linear terms in small parameters apf/er, h/ef 
are taken into account in describing the FSOL. Spon- 
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taneous background spin currents, arising in the two- 
dimensional electron systems with spin-orbit interaction 
without any proximity effects^!, contain higher powers of 
these parameters which are assumed small below. Fur- 
ther, a spontaneous electric surface current, carried by 
Andreev surface states, arises in the superconductor due 
to proximity to the FSOL. Respective structures of wave 
functions and spectra of the surface states are strongly 
influenced by the spin-orbit interaction and the exchange 
field and differ for quasiparticles with opposite momen- 
tum components pf\\ parallel to the surface. The sponta- 
neous electric current arises as a result of this asymmetry 
of Andreev states with respect to Pf\\ — > —Pf\\, to some 
extent analogously to the current induced by chiral sur- 
face or interface states. Conditions for electric and spin 
currents to show up in the system we study, significantly 
differ from each other. Thus, the spontaneous spin cur- 
rent in the FSOL arises even within the framework, when 
the surface electric current vanishes. 

Consider an isotropic s-wave superconductor at x > d, 
deposited with a layer of thickness d made of a ferro- 
magnetic metal. Let a macroscopic thickness of the layer 
be much less than the superconducting coherence length: 
d -C £ s . Both the internal exchange field h and spin-orbit 
Rashba term wcr — a(n x p\\)cr enter the Hamiltonian 

density of the FSOL: H(x) = H (0) - (h(x) + w(x)) cr. 
Here H ^ describes the kinetic energy of free electrons, 
n is the unit vector along the surface normal and p\\ the 
momentum component parallel to the surface. The ex- 
change field is assumed always aligned along the z-axis. 
Both h(x) and w(x) are taken finite and spatially con- 
stant within the FSOL < x < d. The x-axis is taken 
directed into the depth of the superconductor and the 
system is confined by an impenetrable wall at x = 0. 

We assume A <§C h,apf <C £/ and describe the system 
in question by quasiclassical Eilenberger equations for 
Matsubara Green's function: 

do 

-iv fyX — = [(ie n r z + t z A + hs + wf z s) ,g] , (1) 
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g = -7T 



(2) 



Here g(x,pf,e n ) takes 4x4 matrix form in the four- 
dimensional product space of particle-hole and spin vari- 
ables. In the particle-hole space 
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g(p,e n ,x) f(p,e n ,x) 
f(p,e n ,x) g(p,e n ,x) 



(3) 



where all matrix elements are 2x2 matrices in spin space. 
Pauli-matrices in particle- hole space are fj, f± = T x ±iT y , 
while in spin space <Ti. The superconducting order pa- 
rameter matrix is A = l/2[f + A — T^A*]ia y . The opera- 
tor for quasiparticle spin is (l/2)sf z , whereas the opera- 
tor s = l/2[(l + f z )er— (1 — T z )<j y <j<jy\ enters the Zeeman 
term. The order parameter A is taken spatially con- 
stant throughout the superconducting half-space x > d. 



As this follows from recent results for two-dimensional 
superconductors with spin-orbit coupling 30 , a possibility 
for proximity-induced inhomogeneous phase of the order 
parameter in the plane parallel to the interface should 
be studied for sufficiently thin superconducting layer in 
proximity to the FSOL. This two-dimensional inhomoge- 
neous profile of the phase does not appear, however, for 
a massive superconducting sample. 

Electric and spin current densities can be expressed via 
quasiclassical Green's function as follows 



j =N f T{v f Y,SV29{Pf,e n )) i 



(4) 



(•5) 



Here Nf is the normal state density of states per spin 
direction, < ... >g, means averaging over quasiparticle 
states at the Fermi surface. Spin current j s a carries i-th 
spin component along Z-axis in coordinate space. One 
can introduce scalar g Q and vector g components of the 
Green's function in spin space g(pf,e n ) — go{Pf,£ n )&o + 
g(pf,e n )a. As this is seen from Eqs.(@J, (JSJ), go deter- 
mines electric current, while spin current is associated 
with g. 

The Green's function for the FSOL satisfies conven- 
tional boundary conditions on the impenetrable wall at 
x = 0: g(0,pf,e n ) — g(0,pf,e n ), where pf and pf are 
the incoming and the outgoing quasiparticle momenta re- 
spectively. We match solutions of Eilenberger equations 
for the superconduting half-space and for the FSOL with 
the continuity conditions on a transparent interface at 
x = d. Substituting the final result for the Green's func- 
tion into Eqs.QJ, (JSJ, we find no spontaneous electric 
current in the system and finite components j*, j zy of 
spin current situated in the FSOL and flowing parallel to 
the surface. 

One can show that the whole spin current is carried by 
Andreev surface states taking place in the system. We 
find two dispersive branches of Andreev surface states, 
whose energies depend on momentum component parallel 
to the surface: 
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In the absence of spin-orbit interaction spectra of An- 
dreev states, descibed by Eqs. 0, ®, reduce to the 
results for spin-discriminated Andreev states at a ferro- 
magnetic surface^. 

Andreev surface states carry no spin current in a sin- 
glet superconductor, since particles and holes, occupy- 
ing the state, have identical spatially constant local spin 
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polarization and opposite velocities. However, the wave 
function of Andreev surface states does not vanish in 
the FSOL and has a qualitatively different spin struc- 
ture there, as compared with the superconducting re- 
gion. One can extract pole-like terms from the whole 
expression for the electron retarded Green's function 
g R (x,pf,e) near bound state energies £1,2- We deter- 
mine the spin structure of electrons in the Andreev states 
in terms of eigenvectors of these pole-like terms in spin 
space (^)(x,Pf,e). The unit vector P e , describing elec- 
tron spin polarization, can be found from the equation 
P e erQ) (x,pf, e) — As a result, we obtain 

the following spatially dependent spin polarization for 
electrons in Andreev surface states (JHJl at < x < d: 



P e (Pf ,£1,2) 
-(e+ x (e+ 
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(9) 



Local spin polarization of electrons, occupying Andreev 
states, is spatially constant inside the superconductor 
and takes there the same value as follows from Eq.© 
at x = d. Parallel and normal to the surface compo- 
nents of spin polarizations, taken for incoming and out- 
going electrons in one and the same Andreev state, are 
related with each other as P,f(p/, £1,2) = Pf(p/, £1,2), 

p ±(Pf,Ei,2) = -P±(Pf>£l,2)- Also > since £ i = - £ 2, we 
find from Eq.©, that P e (pf,-e m ) = -P e (pf,e m ). 

Spin polarization P h for holes, occupying Andreev 
states, can be derived from Eq.©. The quantity 
P h satisfies the equation —P h ia y &ia y (^)(x,pf,e) — 

(°^)(x,pf,s), which contains the spin operator for holes 

— {l/2)id y aidy. Here (% h )(x,pf,e) is the eigenvector of 

the pole-like term in the Green's function g R (x,pf,e) 
near E\ or £2. As this follows from the general rela- 
tion g R (x,pf,e) — g A (x,—pf,—s), the eigenstates for 
holes and for electrons are associated with each other as 
(^)0,P/, £ ) = (f,)(x,-p f ,-s). Hence, the spin po- 
larization for holes in the state (^ h )(x,pf,e) coincides 

with that for electrons in the state )(x, — p/, — e). 
Further, as this follows from the equation for P e and the 
relation between P e (pf, e m ) and P e (pf, — £ m )> the quan- 
tity P h (pf, e m ) coincides with electron spin polarization 
in the state (_"g„) (x,pf x , — pyy , £,„). Comparing elec- 
tron spin polarizations of states (_"g») (x,Pfx, — Pf\\, £ m) 
and (^j (x,pf x , — Pf\\, £m)i we find finally that spin 
polarization for holes can be found from Eq.© as 
P^P/, £1,2) = P|f(p/x,-P/||,ei,2), Pf(p/,ei,2) = 

-*±(P/s,-P/||,£l,2)< 

The first term in Eq. © describes P£ component of the 
spin polarization, while the second and third terms de- 
termine Pf. Under the transformation pj|| — > —Pf\\ one 



finds e± — > e T and Q± — > T . The first two terms in the 
square brackets in Eq.© are responsible for a spatially 
dependent difference between electron and hole local spin 
polarizations, taking place in Andreev states in FSOL as 
the joint effect of Zeeman and spin-orbit couplings. In- 
deed, for vanishing h or a vectors e± become parallel 
to each other. Then the first and the second terms in 
Eq.© vanish, resulting in identical spin polarizations of 
electrons and holes. Also, at x = d this follows from 
Eq.® P e (p f ,Ei,2)=P h {pf,e 1 ,2). 

Different spin polarizations and almost opposite veloc- 
ities of electrons and holes, occupying Andreev surface 
states Eq.©, result in a net spin current in the FSOL. 
The local spin current density carried by two Andreev 
states can be represented as j?« — + j"'^, where 
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(vf,\\W m [if (p/,£ m ) - if (p/,£ m )] n f {e m )) Sj . 

(10) 

Here W m = ^irANf sin^ is the weight of the delta- 
peak in the local density of states, taken in the FSOL 
for m-th Andreev state (to = 1,2), and n/(e) is the 
Fermi distribution function for quasiparticles. Substitu- 
tion of the represented results into Eq. (|10|) gives exactly 
the spin current density, which follows from Eq.© and 
respective solutions of the Eilenberger equations for the 
Green's function. 

As this follows from Eq. (|10fl after integration over the 
Fermi surface, only parallel to the surface components 
jy z (x) and 2% y (x) of the spin current remain finite. The 
spin current m , carrying along the surface perpendicu- 
lar to the surface spin component, vanishes in accordance 
with the relation P±(fif ,£1,2) — —P±(Pf, £1,2), since the 
contributions from incoming and outgoing electrons, as 
well as holes, cancel each other. The proximity-induced 
background spin current we have found does not lead to 
any spin accumulation and satisfies the continuity equa- 
tion J^i dj s a /dxi — 0. In the limit of small Zeeman cou- 
pling ft < apj, we find the following simple estimations 
for the spin current in the thin layer with spin-orbit in- 
teraction d <C Vf/(apf): 



3ap — ^a/ 



(— 

\ap f 



( dapf 



JcL ■ 



(11) 



Here a, ft = y, z and a ^ 0, A a p > is a constant the 
order of unity, j cL = n s A/p-f is the Landau depairing 
current density. At low temperatures j C L ~ NrVfA. 

In the opposite limit h 3> apf , when the exchange field 
in the FSOL significantly exceeds spin-orbit coupling, es- 
timations for the two components of the spin current give 
different results: 



f yz = B yz ( d ± 

f ap f \ 3 ( dh 



JcL 



(12) 
(13) 



Here B yz , B zy are constants of the order of unity. For 
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h^apf^Vf/d spontaneous spin current densities turn 
out to be the order of j' c j> 

Background spin current density, arising without 
any proximity effects in the two-dimensional metal 
with Rashba spin-orbit interaction^ 7 ., takes the form 
(apf/ef) 3 EfNfVf/6. It is of the third-order in parame- 
ter apf/sf, which is presumably small. These spin cur- 
rents are carried by all the occupied states at a given 
temperatur e) 27 ! 28 , in contrast with the currents induced 
by a proximity to the superconductor. For this reason 
respective reference quantity SfNtVf contains a large pa- 
rameter e/A as compared with j c j,. 

We return now to the problem of spontaneous surface 
electric current. Each separate Andreev surface state, 
taken for given pii, carries finite surface electric current. 
There is no net electric current under the conditions con- 
sidered above, since electric currents carried by Andreev 
surface states Eq.(JBJl, Q, © with p« and — p« cancel 
each other. This is associated with the symmetry of 
scalar component go of the quasiclassical Green's func- 
tion with respect to the sign reversal of the momentum 
parallel to the surface. Spin current takes finite values 
since vector component g of the Green's function does 
not possess the symmetry. However, the symmetry of 
go turns out to be approximate, taking place only un- 
der the conditions apf, h <C £/, within the quasiclassical 
approximation applied to the FSOL. For this reason we 
find below finite spontaneous surface electric current, as- 
suming A <C apf,h < Ef and applying the S- matrix 
approach for describing the FSOL. 

A surface with the FSOL is characterized by the 
normal-state scattering S-matrix, contained reflection 
amplitudes for quasiparticles. The S'-matrix can be 

represented as S = S(l + f z )/2 + 5(1 — f z )/2, where 

S(p n ) = S tr (-p n ) and 



h 



h + w\ 



Here r^i = e lcl 



(14) 

= e ieT -i and, assuming spatially constant h 
and a in the FSOL, 

\Pf* 



©r 



[ 



7T + 2arctan 



■2\Pfx\d, (15) 



where Fermi momenta in the normal metal p f and in the 
FSOL p/t,x satisfy the relation pj x ^ , = p"j x ± 2m\h + 

w (Pf\\)\- 

Making use of explicit expression for the S'-matrix 114|) 
and following the quasiclassical approach with Riccati 
amplitudes in describing the superconductosi^i, we ob- 
tain the quasiclassical Green's function. In particular, we 
obtain spectra of Andreev surface states, which take the 
following form now 



£l )2 = sgn 



sm 



Here X(pf\\) 
®l(-Pf\\)) and 



A=F$ 



A cos 



2 (e T ( P/ ||) + e&fu) - e T (- 

$(p/||) is defind in Eq.©, 



(16) 

■Pf\\) ~ 
where 



one should use the generalized definition for 0±(p/||): 
©±(P/||) = @t( = ' = P/||) — For a small parame- 

ter \h + w\/ef ^ 1 the quantity X(pj\\) vanishes in the 
first approximation, while the definition for Q± reduces 
to that given in Eq.JSJ). 

In general, energies ei^(Pf\\) in Ea. (fC1)|) are situated 
asymmetrically with respect to the Fermi level for a given 
Pfu. Since X(pf\\) and $(p^||) are odd and even func- 
tions of pf\\ respectively, each energy branch £i, 2 (j>/||) 
in Ea. H16|) . as well as the Andreev spectra as a whole, 
is neither odd nor even with respect to the transforma- 
tion p/|| -» -PfW- £i,2(-P/||) = -e 2 ,i(P/||)- As arc- 
suit of the asymmetry, the spontaneous electric current 
density j y (x), flowing along the surface perpendicular to 
the exchange field in the superconducting region, arises 
near the surface with the FSOL. The spontaneous sur- 
face current density at the interface x = d takes compara- 
tively simple form in the case of small spin-orbit coupling 
apj <C (e/ ± h): 



3y{d) = 
_ neNfA 



v fv 



A 

2T 



• 7 ©o , -2 AcosOo 
sin — cosh 



2T 



On , A COS On \ 

ros^tanh 2T ° jX(p n )^ .( 1 1 ) 



Here On, taken for zero spin-orbit coupling, is defined 
as On = 0+| Q= o = 0-U=o- The expression for X{pfit) 
in Ea. l|17[) should be taken linear in small parameters 
apf/(sf±h). Then X{pf\\) cx w z = apf v and averaging 
over the Fermi surface in Eq. I|17|) gives nonzero result for 
j y , while j z vanishes. 

We notice, that an expression for the Josephson criti- 
cal current in S-F-S junctions with small momentum de- 
pendent transparencies D(pf\\)2£ can be obtained from 
Eq.JEJ by replacing X{p n )v fy -» -2D{p n )\v fx \. This 
is not surprising, since both the spontaneous surface cur- 
rent and the Josephson current are actually the two com- 
ponents of total supercurrent carried by the same An- 
dreev interface states, which reduce to surface states in 
the tunneling limit. In the particular case h ~ apf < £j 
the spontaneous surface electric current j y oc apfh/e^ 
is of the second order in a small parameter (h/sf) ~ 
(apf/ef). Since these small second-order terms are dis- 
regarded within the quasiclassical approach to describing 
the FSOL, solutions of Eq.Q) found above show no spon- 
taneous electric surface current, in contrast with the spin 
currents in the FSOL. 
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